Interfacial thermal transport in atomic junctions 
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We study ballistic interfacial thermal transport across atomic junctions. Exact expressions for 
phonon transmission coefficients are derived for thermal transport in one-junction and two-junction 
chains, and verified by numerical calculation based on a nonequilibrium Green's function method. 
For a single-junction case, we find that the phonon transmission coefficient typically decreases 
monotonically with increasing freqency. However, in the range between equal frequency spectrum 
and equal acoustic impedance, it increases first then decreases, which explains why the Kapitza 
resistance calculated from the acoustic mismatch model is far larger than the experimental values 
at low temperatures. The junction thermal conductance reaches a maximum when the interfacial 
coupling equals the harmonic average of the spring constants of the two semi-infinite chains. For 
three-dimensional junctions, in the weak coupling limit, we find that the conductance is proportional 
to the square of the interfacial coupling, while for intermediate coupling strength the conductance 
is approximately proportional to the interfacial coupling strength. For two-junction chains, the 
transmission coefficient oscillates with the frequency due to interference effects. The oscillations 
between the two envelop lines can be understood analytically, thus providing guidelines in designing 
phonon frequency filters. 

PACS numbers: 66.70.-f, 05.60.-k, 44.10. +i, 



I. INTRODUCTION 



In the past decade there has been a significant re- 
search focus on thermal transport in micro scaled. Sev- 
eral conceptual thermal devices, such as thermal recti- 
fiers/diodes, thermal transistors, thermal logical gates, 
and thermal memory^— , have been proposed, which, in 
principle, make it possible to control heat due to phonons 
and process information with phonons. The issue of 
quantum thermal transport in nanostructures was also 
addressed 6 . In this context, the critical information is 
in phonon transmission coefficients that in quasi-one- 
dimensional atomic models can be calculated by trans- 
fer matrix method 7 - - — . However, the evaluation of the 
transfer matrix may be numerically unstable, particu- 
larly when the system size becomes large. Alternatively, 
nonequilibrium Green's function (NEGF) method is an 
efficient way to calculate the transmission coefficient 11 . 
Unfortunately, both of these two methods are numerical 
in nature and do not give analytical expressions. 

For thermal transport and control, the interfacial ther- 
mal scattering process is becoming increasingly impor- 
tant, especially in practical devices. Two theories, acous- 
tic mismatch model^ and the diffuse mismatch modet^, 
have been proposed to study the mechanism of the ther- 
mal interfacial resistance. However, both models of- 
fer limited accuracy in nanoscale interfacial resistance 
predictions^ because they neglect atomic details of ac- 
tual interfaces. A scattering boundary method within the 
lattice dynamic approach was first proposed by Lump- 
kin and Saslow to study the Kapitza conductance in 



a one-dimensional (ID) lattice^, and was then applied 
to calculate the Kapitza resistance in two- and three- 
dimensional (3D) lattice s 16 ! 17 . This method can predict 
thermal interfacial conductance between heterogeneous 
materials with full consideration of the atomic struc- 
tures in the interface. Recently, this method was ap- 
plied to study the ballistic thermal transport in nan- 
otube junctions 18 , spin chains^, and honeycomb lattice 
ribbons^. 

In this paper we give an explicit analytical expression 
of transmission coefficient obtained through the scatter- 
ing boundary method, and use it to study the interfa- 
cial thermal transport across atomic junctions. First, in 
Sec. HU we introduce a model in which two semi-infinite 
ID atomic chains are coupled either via a point junction 
or an extended junction region. By using the bound- 
ary scattering method we derive the exact expressions 
for phonon transmission coefficients for thermal trans- 
port in one-junction and two-junction chains in Sec. IIIII 
The role of various parameters on the junction conduc- 
tance is analyzed and discussed in Sec. [IVJ In section 
IIVI we also estimate the interfacial conductance between 
two 3D solids. In Sec. E] we introduce briefly the NEGF 
method, and use it to verify the results from analytical 
formulae for the thermal transport in our model. A short 
summary is presented in Sec. IVI1 



II. MODEL 

The one-dimensional atomic chain consists of three 
parts: two semi-infinite leads and an center region (see 
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FIG. 1: (color online) A schematic representation of the ID 
atomic chain model. The size of the center part is Nc = 
8. The left and right regions are two semi-infinite harmonic 
atomic chains at different temperatures Tl and Tr. The three 
parts are coupled by harmonic springs with constant strength 
ki2 and £23; all of which are harmonic chains with mass and 
spring constant as mi, ki, mi, k-z and 7713, /C3, respectively. 



Fig. Q}. The two leads are in equilibrium at different 
temperatures Tl and Tr. The three parts are coupled 
by harmonic springs with constant strength k±2 and ^23; 
all of which are harmonic chains with mass and spring 
constants mi, ki, TO2, fc 2 and m.3, fca, respectively. So the 
total Hamiltonian can be written as 

H= ^ H a + ^k 1 2(x 1A -X2.i) 2 + ^k 2 3{x2,N t! -X3,i)' 2 ; 



q=1,2,3 



here, 



1=1 



N a -1 



(i) 



) 2 . (2) 



Where x a j is the relative displacement of i-th atom in 
a-th part. If there is no center part, that is, the two 
semi- infinite leads connected directly by ki2, then by set- 
ting a = 1,2 and fc 2 3 = in Eq. (Q}, we can obtain the 
corresponding Hamiltonian. For the semi-infinite leads, 
N a = 00. 



We first consider a point-junction case, that is, two 
semi-infinite harmonic chains connected by a spring with 
constant strength ki2- We assume a wave solution trans- 
mitting from the left lead to the right lead. We la- 
bel the atoms as —00, • • • , —1, 0, 1, 2, • • • , +00. Atoms 
and 1 are connected by ki2 spring. An incident wave 
from left is assumed as xi — X{e~~ lult . When it ar- 
rives at the interface, it will be partially reflected and 
partially transmitted. The reflected wave amplitude is 
xr = ri2\i J e~ lut and the transmission wave can be 
written as xt = t\2\^ 2 1 e 

(Af 1 + r l2 Xi)e- iut , x = (1 



£i 2 A 2 ±e tuJt - So at each atom we have 



X-l 



r\2)e~ 



xi = t 12 e-' lult , x 2 = ii 2 A 2 e- laJt , Here, Xj = e lq ? a \ 
qj is the wave vector, aj is the interatomic spacing. For 
the atom in the j — th part, we can have the equation of 
motion as 



drx-j 



dt 2 



kjixj^n+l Xj^n) ~f" ^j{xj : n Xj^n — l), (5) 



each wave transport separately and satisfies such equa- 
tion. Thus Aj satisfies the dispersion relation of the cor- 
responding lead as 



uj 2 mn 



- kj Xf 



2k, 



kjXj. 



(6) 



The quadratic equation has two roots. Which one should 
we choose? Replacing uj with uj + irj, rj = + , none of the 
eigenvalues A will have modulus exactly 1. We find for 
the traveling waves^i 



= !-»?-, 

v 



(7) 



thus the forward moving waves with group velocity v > 
have |A| < 1. Therefore we should take the one with 
|A| < 1 of the two roots which are given as 



III. ANALYTICAL SOLUTION FROM THE 
SCATTERING BOUNDARY METHOD 



-hi± 



h J = -^(uj + ir 1 ) 2 -2. (8) 



Heat current flowing from left to right through a junc- 
tion connecting two leads kept at different equilibrium 
heat-bath temperatures Tl and Tr is given by the Lan- 
dauer formula 6 
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hu: - f R (uj)]T[uj}duj, 



(3) 



which allows us to develop the junction conductance for- 
mula 
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2~tt 



— — / du huj T \uj 



OT ' 



(4) 



here, Jl,r = {exp^w/^BT^)] — is the Bose- 

Einstcin distribution for phonons, and T[uS\ is the fre- 
quency dependent transmission coefficient. Therefore, 
the key step for the thermal transport characterization is 
to calculate the transmission coefficients. 



From the scattering boundary method, the coefficients 
T\2, ti2 can be obtained from the continuity condition at 
the interface as: 



u} 2 m 1 x a = —kiX-i + (fci + fci 2 ).x - ki 2 xi; 
uj 2 m 2 xi = -ki 2 x a + (ki2 + k 2 )xi - fci 2 x 2 . 

Finally we can get the transmission coefficient as 
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(9) 
(10) 



(11) 



here. 



h{Xi - l/\i)(kj 



kj/Xj) 



{ki kij ki j Xi)(kj 



(12) 

Of course, we can also use t\2 to express T[u] as 
%^\t 12 \ 2 , here the group velocity v t = % = 



3 



'4k, 



uj 2 , which is derived from the dispersion rela- 
tion given by Eq. Thus, the transmission coefficient 
can also be expressed as 



TM 



2m 2 



mi 



\\tl 



(13) 



here 



— kijki(Xi — 1/Aj) 



kij hi/ \i)\kj kij kj/Xj) 



k 2 ' 



(14) 



For the long- wave limit, that is, ui — + , we get r 

\fk~ii 



^/kimi + ^J kjirij 1 



i J 



and the transmission is 



T[0 + ] = 



Ay / kimik 2 m2 



{y/kimi + \fh2jn2) 1 



(15) 



This result is consistent with the one obtained for the 
acoustic mismatch model, i.e., T — rg^ffp - 2=2- Where 
the acoustic impedance is Z% = piVi = (mi/ai)i)i, and 
Zi(u) = + ) = yjkirrii. We note that in acoustic mis- 
match model the transmission coefficient is frequency in- 
dependent, and in reality it only applies in the limit of low 
frequency/long wavelengths. In this case the phonon sees 
the interface only as a discontinuity between two semi- 
infinite media and the transmission does not depend on 
the coupling spring strength fey. If the two leads have 
the same acoustic impedance for long wave limit, then 
T[0+] = 1; otherwise T[0+] < 1. 

For a two-junction case, which is shown in Fig. [TJ the 
transmission wave will be reflected and transmitted by 
the second boundary, leading to multiple reflections. Fi- 
nally the total transmitted wave function is obtained as 
a superposition of multiple reflections and transmissions, 
resulting in the transmission coefficient through the cen- 
ter part 



TM 



(l-|r 12 | 2 )(l-|r 23 | 2 ) 

h ,2(JV C -1)|2 : 

|1 - 7-237-21A2 I 



(16) 



here ry and A.; are determined by Eq. (fT2")l and Eq. (|5J|; 
Nc is the number of atoms in the center atomic chain. 
From this expression, we can find that the transmis- 
sion coefficient oscillates with frequency, and is between 
the envelope lines of maximum and minimum transmis- 
sion, which are T m M = (1 - |r 12 | 2 )(l - |r 23 | 2 )/(l - 

interference and T m i n [cj] = 



i-maxM — 

1^23^21 1) 2 for constructive 



/(i-ki 2 | 2 )(i 

terference. 



1^23 1 2 ) / (1 + |^23^2i|) for destructive in- 



IV. RESULTS AND DISCUSSIONS 



A. Thermal transport in ID one-junction chains 




FIG. 2: (color online) The transmission coefficient vs fre- 
quency uj for different interface coupling k\2 in one-junction 
chains. (a) shows the transmission in one junction con- 
nected by the same semi-infinite atomic chains with ki = 
&2 = 1.0, mi = m.2 = 1.0; the solid, dashed, dotted and 
dash-dotted lines correspond to k\2 = 0.1, 0.5, 1.0 and 2.0, 
respectively. (b) shows the transmission in one junction 
connected by two different semi-infinite atomic chains with 
ki — 1.0, mi — 1.0, hi = 3.0 and 7712 = 4.0; the solid, 
dashed, dotted, dash-dotted and shot-dashed lines correspond 
to fei2 = 0.5, 1.0, 1.5, 3.0 and 8.0, respectively. 



cases Eq. (fTTj) . Eq. (|T2")) and Eq. (fTBj) by using the scat- 
tering boundary method. Using these analytical expres- 
sions, we analyze the role of various parameters on the 
thermal transport in one- and two- point junctions. 

Figure [2] shows the transmission coefficient as a func- 
tion of frequency for a different interface spring constant 
feia for the point-junction model. The maximum fre- 
quency at which the transmission coefficient is above zero 
is equal to the minimum of 2^/ki/mi and 2y / fc 2 /m 2 . 
In Fig. [21a), the two semi- infinite atomic chains have 
the same mass and spring constant. When the inter- 
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In Sec. IIIIl we have derived the analytical expres- 
sions for the phonon transmission coefficient for point- 
junction and extended-junction (two point junction) 



FIG. 3: (color online) The thermal conductance vs interface 
coupling fci2 in point-junction model. Here, ki = 1.0, mi = 
1.0. 
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FIG. 4: (color online) The thermal conductance vs the ra- 
tio of fci2/fci2m in one-junction atomic chain. Here fci2 m is 
the harmonic average of the spring constants of the two semi- 
infinite leads, (a) k\ = 1.0, mi = m% = 1.0; the solid, dashed, 
and dotted lines correspond to = 0.1, 1.0, and 40.0, respec- 
tively, (b) fci = 1.0, mi = 1.0, ki = 10.0; the solid, dashed, 
and dotted lines correspond to m2 = 0.01, 1.0, and 100.0, 
respectively. 



face coupling k\2 equals to that of the chains, the trans- 
mission is equal to one in the whole frequency domain, 
because of the homogeneity of the chain structure. If 
fci2 increases or decreases, the transmission coefficient 
decreases. If we set k\jm\ — kijmi, the transmission 
coefficient exhibits similar behavior, the only difference 
is that the transmission coefficient changes to the value 
obtained by Eq. (fT5|i . In Fig. [^fb), the two semi-infinite 
atomic chains have different masses and spring constants. 
The transmission decreases with increased frequency for 
all the coupling values k\2- Also, it appears that for a 
given frequency the transmission is maximized for a k\i 
value residing between k\ and k2- From Eq. (jlip and 
Eq. ([12]), T[u] = 0, if k u = 0; and T[u>] has definite 



value 1 



1 ktii-x-^+^ii-x- 1 ) 1 



if k w = 



The maximum transmission concept results in the 
maximum junction conductance as shown in Fig. [21 With 
the increasing of fci2, we find that the conductance will 
first increase, then arrive at maximum value, and then 
slightly decrease and at last it will tend to a constant. 
We find that the maximum transmission or conductance 
occurs at ki2 given by 



ki2 = k 12m = 



2fcifc 2 
ki_+k 



(17) 



i.e., when the coupling spring stiffness is equal to the 
harmonic average of spring connecting atoms in the two 
semi-infinite chains. In Fig. 21 we show the thermal con- 
ductance vs the ratio of fci2 and k\2 m . For the two semi- 
infinite chains with the same mass m\ = m,2, the max- 
imum conductance occurs exactly at k\2m- If the two 
leads have different masses m\ ^ rri2, the maximum con- 



FIG. 5: (color online) The transmission coefficient vs fre- 
quency for different mass ratios mi/ mi at the interface cou- 
pling fci2m- Here, ki = 1.0, ki = 3.0, fci2 = fci2m = 1-5 and 
mi = 1.0. 



ductance is almost exactly at the fci2 m point, for mass 
ratios ranging from 0.01 to 100. 

In Fig. [5l we show the curves of the transmission as 
a function of frequency for interface coupling equal to 
k\2m- If fci/mi = k2/ni2, that is, when both chains have 
the same frequency spectrum of [0, 2^k\/m\], the trans- 
mission equals to a constant T[uj] = T[0 + ], which can 
be seen from the solid line in Fig. [SJ and which is consis- 
tent with Fig.[2fa). Thus for chains with matched spectra 
the transmission is frequency independent. Let us now fix 
ki, k2 and k%, and decrease rri2- In the range between the 
point of equal-spectrum {u m = k\jm\ = ^2/^2) and the 
one of equal-impedance (Z(u> = + ) = k\m\ = /C2TO2), 
the transmission will first increase with frequency and 
then decrease. Otherwise, there is a monotonic de- 
crease. The former behavior is quite interesting, as one 
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FIG. 6: (color online) The transmission coefficient vs fre- 
quency for different interface coupling fci2m- Here, fci = 
1.0, mi = 1.0, k 2 = 0.7, m 2 = 0.3. 
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FIG. 7: (color online) (a) The cutoff frequency vs interface 
coupling for ID one-junction atomic chains. The parameters 
are: ki = 1.0, mi = 1.0. (b) The transmission as function 
of interface coupling for ID one-junction atomic chains. The 
parameters are: k\ — 1.0, mi = 1.0, ki = 0.7, = 0.3 



expects that the transmission should be the largest in 
the long wavelength limit. For highly dissimilar materi- 
als, the transmission coefficient in the whole frequency 
range is much larger than that in the long wave limit 
T[uj — + ] = (z^+zl) 2 > * nus ^ ne rea, l conductance is far 
larger than that calculated from the acoustic mismatch 
model. This result explain why the interfacial resistance 
calculated from the acoustic mismatch model is far lager 
than the experimental value measured at low tempera- 
tures, where the phonon transport can be regarded as 
ballistic transport. 

In many real interfaces, interface coupling is very weak, 
that is, the k 12 is less than ki 2m . So it is desirable 
to study the thermal transport in atomic chains in the 
weak coupling limit. Figure |H] shows the transmission co- 
efficient as function of interface coupling. In the weak 
coupling limit, with the frequency increasing, the trans- 
mission decreases rapidly to zero, so the frequency region 
where phonons are effectively transmitted is very narrow. 
With interface strength increasing, more and more modes 
contribute to the transmission and the phonon transmis- 
sion window widens. If the interface coupling increases 
further, that is fci 2 /fci 2m > 0.1, out of the weak inter- 
face coupling limit, all the phonons contribute to the 
transmission. The only further change with increasing 
ki2 is the actual values of the transmission coefficients 
increase. In Fig. [Zta), we show the transmission cutoff 
frequency as function of the interface coupling. Here, 
we define the cutoff frequency uj cuto g at which the trans- 
mission r(w cuto ff) = 0.1T(0+). We find that the cutoff 
frequency shows linear dependance on interface coupling 
in the weak coupling limit fci2 < 0.1fci2 m - If the interface 
strength increase further, the cutoff frequency is satu- 
rated. In Fig. EJb), we show the transmission as func- 
tion of interface coupling for several different phonons. 
We find that in the weak interface coupling region, the 
transmission is proportional to the square of the interface 
coupling, which is consistent with the formulas Eq. (|13[) 
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FIG. 8: (color online) The thermal conductance vs interface 
coupling for ID point-junction atomic chains. The parameters 
are: k\ = 1.0, mi = 1.0. 
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FIG. 9: (color online) The thermal conductance vs interface 
coupling for 3D one-junction atomic chains. The parameters 
are the same with Fig. [8] (a) Interface coupling is far less 
than the coupling fci2 m : fci2m/fci2 = 0.001 — 0.1; (b) Interface 
coupling is in the region of 0.1fci2m ~ 0.9fci2m- 



and Eq. (|14[) . In the weak interface coupling region, for 
the ID atomic one-junction chains, it is shown that the 
thermal conductance is linear with the interface coupling 
(see Fig. If we strengthen the interface coupling be- 
tween the two chains, the conductance will be linearly 
enhanced. For different mismatched chains, the absolute 
values of the conductance arc different, but dependence 
on the coupling strength is the same. 



B. Thermal transport in 3D single-interface 
structures 



The thermal conductance Eq. ((4]) can also be written 

.,23. 



p oo 

a = I dui huj T\uj\ 
Jo 



g/H 

dT 



v(uj)D{uj), (18) 



because of v(u>) = duj/dk and phonon density of states 
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FIG. 10: (color online) The transmission coefficient of the 
two-junction atomic chains. Parameters: ki = 1.0, mi 
I J). k 2 = 0.9, m 2 = 1.6, k 3 = 4.5, m 3 = 2.0, The solid, 
dotted, dashed and shot dashed lines correspond to max- 
imum transmission, minimum transmission, N c — 4 and 
N c — 9, respectively. The interface couplings are different: 
(a) fcia = 0.3, k 23 = 0.7; (b) k 12 = 1.0, k 23 = 4.5. 



in ID structure, D(u) = l/(2irv), we can obtain Eq. (j4j. 
In order to estimate the behavior of the interfacial ther- 
mal transport across interfaces in 3D structures, we only- 
need to change the phonon density of states in the above 
equation. Because the density of states for 3D structure 
within the Debye approximation is D(u>) ~ cj 2 , therefore 
we can replace w with w 3 in Eq. ([5]); the thermal con- 
ductance as a function of the coupling strength is shown 
in Fig. [9] From Fig. [9^a) , we find that in the weak in- 
terface limit, conductance is proportional to the square 
of interface coupling, which is consistent with the results 
from other models 2 - - — , while it is linear dependent on 
the interface coupling in ID junctions. This is due to 
the fact that in 3D low frequency region contributes rel- 
atively little to the conductance as the density of states 
is low there. If the interface coupling increases further, 
that is ki2/ki2 m > 0.1, out of the weak interface coupling 
limit, all the modes contribute to the transmittance, the 
conductance is no longer proportional to the square of 
the interface coupling, and the slope continuously de- 
creases. In some intermediate ranges the conductance 
is approximately proportional to the interfacial coupling 
(see Fig. EJb)), which is consistent with the results from 
molecular simulation approach 27 . For stronger coupling 
the conductances for the ID case and 3D one have similar 
behaviors, the slope of both cases will decrease continu- 
ously to be zero at point ki2 m , where the conductance 
will be maximized and then decrease slightly to a limit- 
ing value. 



C. Thermal transport in extended junctions 

Now we focus on a case where the junction is extended 
and involves a center part. The overall behavior of the 




FIG. 11: (color online) The transmission coefficient of the 
two-junction atomic chains. Here, ki = 1.0, mi = 1.0. 
The solid, dotted, dashed and shot dashed lines correspond 
to maximum transmission, minimum transmission, N c — 4 
and N c = 9, respectively, (a) k 2 = 3.0, m 2 — 5.0, k 3 = 
1.0, m 3 = 1.0, fei2 = k 23 = 1.0; (b) k 2 = 3.0, m 2 = 1.0, k 3 = 
5.0, m-3 = 1.0, fci2 = ki 2m = 1.5, k 23 = k 23m = 3.75; (c) 
k 2 = 3.0, m 2 = 3.0, k 3 — 5.0, m 3 — 5.0, ki 2 — k 12m = 
1.5, k 23 = k 23m = 3.75. 



transmission is the combination of the transmission be- 
havior in single point-junction case and the oscillatory 
behavior due to phonon interferences arising form multi- 
ple scattering. We show the transmission coefficient as a 
function of frequency of an arbitrary case in Fig. UOl a) . 
Here, the three chain parts have different masses and 
spring constants, and the interface coupling is not spe- 
cial. From the analytical expression of Eq. (|16l) , we plot 
curves of the maximum transmission and minimum trans- 
mission, N c — 4 and iV c = 9. The transmission oscil- 
lates between the envelop lines of maximum and mini- 
mum transmission. The maximum transmission line will 
increase first, and the minimum transmission line will 
monotonically decrease with frequency. However for in- 
terface coupling that is the same with the leads, the two 
envelop lines will monotonically decrease, which can be 
seen in Fig. flOTb). 

For some special cases, the transmission coefficient in 
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FIG. 12: (color online) The maximum and minimum trans- 
mission coefficient of the two-junction atomic chains. Here, 
ki = 1.0, mi = 1.0. The solid, dashed lines correspond 
to maximum transmission and minimum transmission k 3 
1 .0. m,3 = 1-0, respectively; the dotted and dash-dotted lines 
correspond to maximum transmission and minimum trans- 
mission k 3 = 5.0, 7713 = 5.0, respectively. The inset shows 
the transmission coefficient with frequency for different fe. 
k 1 — k-i = 1.0, m 1 = 7713 = 1.0. The dotted, dashed, and 
solid lines correspond to k 2 = 0.5, 0.1, and 0.02 respectively. 
For all the curves, 7712 = k 2 and ki2 = fci2m, k 23 = k 23m . 
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FIG. 13: (color online) The comparison of the results from 
scattering boundary method and nonequilibrium Green's 
function method for the transmission coefficient in two- 
junction atomic chains. The square curve and solid line cor- 
respond the parameters: N c = 6, ki = 1.0, mi = 1.0, k 2 — 
1.5, m 2 = 1.3, k 3 = 2.0, m 3 = 1.7, fei 2 = 1.3, k 23 = 0.8; 
the circle curve and dashed line correspond the parameters: 
N c — 13, ki = 1.0, mi = 1.0, k 2 = 1.5, m2 = 1.3, k 3 = 
4.0, 777.3 = 2.7, fci2 = 1.3, k 23 = 0.8. The square and circle 
curves are the results from nonequilibrium Green's function 
method; The solid and dash lines are the results from scat- 
tering boundary method. 



the frequency domain has interesting phenomena, which 
are shown in Fig. 111! In Fig. Illf a). the transmission for 
the case of two identical leads is shown. In this case, the 
maximum transmission is equal to one, the infinite-long 
wavelength phonon and the resonance mode can transmit 
fully through the center part. The minimum transmis- 
sion is very low, indicating efficient destructive interfer- 
ence. Figure ITTT b) shows the transmission when all three 
parts are different and connected by interface couplings 
k\2m and fc23m- We find that overall trend for the maxi- 
mum and minimum transmission lines is increasing first, 
then decreasing. If, in addition, the ratios of ki/rrii are 
the same for three parts, then the maximum and mini- 
mum transmission are constants in the whole frequency 
range, and the transmission coefficient through finite-size 
center part oscillate between the two constants, which 
can be clearly seen in Fig. ITTTc). Therefore, we can use 
the above properties of transmission to design the fre- 
quency filters. Figure [12] shows the maximum and mini- 
mum transmission coefficient for the filter. If the spring 
constant of the center part is very different from the ones 
of the the two leads, the oscillatory peak is sharp, and 
transmission for most of the frequency will tend to zero, 
only few resonant frequency can be transmitted. This 
finding provides guidelines for the design of selective fre- 
quency filters. 



V. VERIFICATION BY NONEQUILIBRIUM 
GREEN'S FUNCTION METHOD 

The NEGF method is an exact approach to study the 
ballistic thermal transport through junctions. Following 
the discussion in Sec. UH if we use a transformation for 
the coordinates, Uj — ^/fnjxj, which is called the mass- 
normalized displacement, then the Hamiltonian can be 
written as 

H = ]T H a + UjV 0>2 U 2 , (19) 

a=l,2,3 0=1,3 

where H a = i (Pj P a + £7 J K a U a ) . K a is the mass- 
normalized spring constant matrix, and V12 = (V2i) T is 
the coupling matrix of the left lead to the central re- 
gion and similarly for V23 is the coupling matrix of the 
right lead to the central region. As stated in Ref. ? , 
the element of the coupling matrix V a g is equal to 
—kij/sqrtmirrij which corresponding to the coupling be- 
tween the ith atom in region a and the jth atom in region 
P- 

We can use the nonequilibrium Green's function 
method^ to study the thermal transport in the atomic 
chain. We define the contour-ordered Green's function 
as 

G^{Ty)^- l -(TU a {T)U p (T') T ), (20) 

where a and /? refer to the region that the coordinates 
belong to and T is the contour-ordering operator. Then 
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the equations of motion of the Green's function can be 
derived. In particular, the retarded Green's function for 
the central region in frequency domain is 



G r [uj] = (lo + i V ) 2 -X 2 -E>] 



(21) 



Here, = an d S Q = V2, a g a Va,2 is the self- 

a=l,3 

energy due to interaction with the heat bath, g r a = [{uj + 
ir\) 2 — Ka]' 1 . And in the advanced Green's function G a — 
(G r y , the transmission coefficient can be calculated by 
the so-called Caroli formula as 



Tp a [u]=Tr(G r T p G a T a ) 



(22) 



where r a = »(££[w] - [w]). 

For single-junction atomic chains, if we regard the two 
atoms in the interface (atom and atom 1) as the center 
part, then we can still use the formulae above to study 
the phonon transmission leading to the exact formula 
yielding the same result with the one obtained from the 
scattering boundary method. In Appendix A, We give 
the analytical proof of this fact. 

For two-junction atomic chains, according to the 
NEGF formulas, we do the numerical calculation and plot 
the curves of the transmission coefficient as a function 
of frequency and compare them to the results obtained 
the scattering boundary method (see Fig. 13). We find 
that for any arbitrary case, the results from the NEGF 
method and the scattering boundary method are exactly 
the same. If there is no many-body interaction, that 
is, for the ballistic thermal transport the scattering ma- 
trix approach and the Green's function method give the 
same results. These two methods are equivalent, which 
has been proved from other points of view in Refs i 28 i 29 . 



constants of the leads. At the point near k±2 = fci2m> 
the transmission T[uj] is a constant if kijmi — k\jm\\ 
if not equal, in the range between ki/mi = A:2/?Ti2 and 
k\wi\ = fc 2 TO2, the transmission coefficient increases first 
then decreases with the increasing of frequency, other- 
wise the transmission monotonically decreases as the fre- 
quency increasing. For weak interface coupling, the cutoff 
frequency and the interface conductance for ID chain is 
linear dependent with the interface coupling strength. 

Because of different density of states, we change the 
formula of conductance to mimic the thermal transport 
in 3D junctions. In weak interface coupling limit, we 
find that the conductance is proportional to the square 
of the interface coupling, which is consistent with the re- 
sults from other models. The slope of the conductance 
as function of interfacial coupling strength decreases con- 
tinuously from two to zero, in certain range of which, 
the conductance is linear proportional to the interface 
coupling, which are consistent with the results of other 
molecular simulations. 

For two-junction case, the transmission will oscillate 
with frequency in the envelop lines of maximum and 
minimum transmission which are determined by the one- 
junction picture. The transmission sometimes oscillates 
between two decreasing envelop lines, sometimes between 
two increasing envelop curves, or between two constants, 
etc. 
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VI. CONCLUSION 

In this paper, we study the ballistic interfacial thermal 
transport in atomic junctions, we give the analytical sim- 
ple formulae Eq. JITJ), Eq. JI3J) and Eq. (JTHJ) for the trans- 
mission of one-junction and two-junction cases, which are 
consistent with the results from the NEGF method. 

For one-junction case, we find the transmission and 
conductance are maximized when the interface spring 
constant equals to the harmonic average of the two spring 



Appendix A: Analytical proof of the equality of the 
two methods for one junction 

In this appendix we give the analytical proof for the 
equality of the scattering boundary method and the 
non-equilibrium Green's function approach for the one- 
junction atomic chains. 

From the scattering boundary method, we obtain the 
transmission Eq. (IT31 and Eq. (TTJ} , that is 



T[u] 



\J 4fc2?Ti2 — ui 2 m 



-fci 2 fei(Ai - 1/Ai) 



\J Ak\m\ ~ u 2 m\ 1 (h - ki2 - k 1 /X 1 )(k 2 - k 12 - fc 2 /A 2 ) 



k 2 

ft 12 



(Al) 



From the dispersion relation Eq. ([5]), we can obtain ; and 

kj-kj/X j = u} 2 m j -k j (l-X j ) (A2) k]\Xj - 1/Aj| 2 = ^(ikjrrij - w 2 m|) (A3) 
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, So we can get 



T[u>] 



k\ 2 uj 2 ij 4fcirni — ui 2 m 2 \J Ak 2 m 2 — Lu 2 m 2 
[w 2 mi - fci(l - Ai) - fci2][w 2 m 2 - k 2 {\ - A 2 ) - ku] - kf 2 



(A4) 



Using the NEGF formulae, we regard the two atoms in 
the interface (atom and atom 1) as the center part 0, 
then the dynamic matrix of the center as 



K 



fcl+fcl2 

mi 
-fcl2 



-fci2 

fcl2 + fc2 



(A5) 



And the coupling matrices between the leads (parts 1 
and 2) and the center (part 0) are Voi = (fci/m-i , 0) T 
and V02 — (0 , k 2 /m 2 ) T , and according to Ref^S, we can 
obtain the surface Green's function as 



fJ, 



(A6) 



here, i = 1,2 corresponds to the left and right lead. Then 
we can get the self energy (E r = Vbi<?[Vio + ^025^20) as 



S r = 




(A7) 



Thus we can calculate the retarded Green's function of 
the center G r = (oj 2 I — Kq — which reads as 



G r = 



Ai B 

B An 



j_ f A 2 -B 
A\-B A x 



(A8) 



here, I is two-dimensional identity matrix and 



A 2 k * , x fc 12 

Ai = lo (1 - Aj) : 



B = 



y/rnim 2 



A = A X A 2 - B 2 . 



(A9) 
(A10) 



The advanced Green's function G a equals to (G r )^. And 
from the self energy we can get 



Ci 



Tl I j ' r2 ~ V C 2 







(All) 



here, Ci = ^- yj 'Akiirii — us 2 m 2 . Therefore, we can calcu- 
late the transmission coefficient from the Caroli formula 
Eq. (|22|) . at last we obtain 

B 2 C\C2 B 2 C\Ci 



T[oj] = Tr(G r r 1 G a r 2 ) = 



AA* \A X A 2 -B 2 \ 2 

(A12) 

Inserting the values of Ai,B and Ci, we get exactly the 
same result with Eq. (|A4I) . Therefore, the results from 
the scattering boundary method and non-equilibrium 
Green's function approach are equivalent. 
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